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Abstract 

The transmission of two electrons through a region where they interact is 

found to be enhanced by a renormalization of the repulsive interaction. For a 

specific example of the single-particle Hamiltonian, which includes a strongly 

attractive potential, the renormalized interaction becomes attractive, and the 

transmission has a pronounced maximum as function of the depth of the 

single-electron attractive potential. The results apply directly to a simple 

model of scattering of two interacting electrons by a quantum dot. 
73.23.-b, 72.10.-d 
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Transport through quantum dots, |J systems of quantum dots, @ or one-dimensional 
wires [[J is currently the subject of much interest. Such structures have potential applica- 
tions as artificial atomic or molecular devices, and are very instrumental in studying strong 
correlations and their implications, e.g., the Kondo effect 0] or magnetic transitions. || 
The theoretical studies of the transmission through such systems concentrate mostly on the 
single-electron conductance, |J with the electron-electron interaction taken into account in 
an approximate manner, e.g., via the charging energy (introducing the notion of the dot 
capacitance [0]) or by using the slave-boson mean-field approximation. || Exact solutions 
of interacting electrons, in this context, focus mainly on the electronic spectrum of two 
electrons, confined to a restricted region in space, which models the quantum dot (see, e.g., 
Ref. ||). Obviously, two is the minimal number of electrons required to study interac- 
tions. However, devices including two electrons can be investigated experimentally, flJlTi 



and may have potential applications in quantum computers. flU]] Recently, we have studied 



the electronic spectrum of two interacting electrons on artifical atoms, [EJ and found that 
the understanding of the full physical behavior of such systems requires a consideration of 
the entire system, including the leads. In that calculation we found that interactions can 
delocalize one of the electrons, depending on the strength of the connections between the 
dot and the leads. In other words, it is not sufficient to treat the interacting region as an 
isolated system. This transition of an electron from the bound state into the band may be 
relevant for the metal-insulator transition brought about by donor ionization. The exact 
solution of the two electron problem yields a renormalized value of the interaction, which 
depends on the energy. 

In this paper we study the one-dimensional transmission of two electrons, when they 
interact via a contact potential at a certain point. We consider the Hamiltonian 

H(xx, ax, x 2 , a 2 ) = TL {xi,x 2 ) + U8{x 1 )8{x 2 )8 (Tl - a2 , 

Ho(x 1 ,x 2 ) =H sp (x 1 ) +H sp (x 2 ), (1) 

where Xi and cr, are the coordinate and the spin component of the ith electron, and 7i sp is 



the single-particle Hamiltonian, independent of the spin components. In ([!]), U is the local 
(Hubbard) interaction. This is a continuum version of the Anderson impurity model. [13| 
Since 7io does not depend on the spins, and since the interaction vanishes on the triplet 
states (which are antisymmetric in space and vanish when x\ = 22 = 0), it is convenient 
to separate the Hilbert space into triplet and singlet states. The former are unaffected by 
the interaction, and their scattering is fully described by the single-electron non-interacting 
transmission. The rest of this paper is devoted to an exact solution of the singlet case. We 
express the scattered wave function in terms of the single-particle spectrum, and calculate 
the two-electron current. 

Our explicit calculations are carried out for a single-particle Hamiltonian which has an 
attractive short-range potential, resulting with a bound state, whose inverse localization 
length k determines the electric current (or equivalently the transmission) in the absence 
of the interaction. When one views the model as a 'quantum dot' coupled to ideal one- 
dimensional conductors ('leads'), then the strength of the attractive single-particle potential 
can be regarded as the gate voltage applied to the 'dot'. |I],[2| 

We consider explicitly two scenarios for the transmitted current. In the first case, the 
incident wave contains two propagating electrons (with wave numbers p\ and p 2 ). For our 
contact interaction [Eq. (1)], we show that in this case there is no effect of the interaction: 
for a sample of large length L, the current is given by 

J = e( Pl \t pl \ 2 +p 2 \t P2 \ 2 )/(mL) + 0(U/L 2 ), (2) 

where t p is the single-electron non-interacting transmission amplitude. The reason for this 
simple behavior is quite clear: In our model, the two electrons "feel" one another only when 
both are at x = 0, and the amplitude for this to occur is proportional to 1/L. 

In the second case, the "impinging" wave contains one propagating electron (with wave 
number p), while the other is captured by the attractive potential. Now the current is found 
to have the form 

j = epT(U,p)/(mL), (3) 



and the effective transmission T(U,p) has a very interesting dependence on the single- 
electron parameters and on U. Specifically, Fig. [I] shows the dependence of T(U,p) on k, 
for several values of U and p (in units of h 2 /{2m) and of the large wave number cutoff u, 
respectively). Clearly, T(U,p) increases significantly when U > (compared to the case 
U = 0), reflecting the "screening" of the attractive single-electron potential by the bound 
electron, eventually leading to a transition from the doubly bound 'insulating' state to a 



state in which one electron is 'free', which can be called 'metallic'. |12fl For the special 
S— function attractive potential discussed below, and for large k, the arguments of Ref. [l2j 
imply that this transition happens for U > 2. Indeed, when U < 2 Fig. [I] (left panels) 
shows a monotonic decrease in T, but T is much larger than the non-interacting result. At 
U = 2 (central panels), T is found to approach a finite plateau value even at large k. This 
value increases with p. For U > 2, we observe a peak in T, which moves to larger values 
of k as U approaches 2 from above. This peak has T»l for small p, but it decreases and 
broadens as p increases. As we show below, this peak corresponds to the resonance with 
the just-bound doubly-occupied state. T(U,p) decreases for larger k, when the effective 
renormalized interaction becomes more and more attractive. 

In addition to the above peak (and to the peak at k = 0), we sometimes find a third 
peak at p ~ k <C u, i. e. when the total energy of the two electrons vanishes; at larger k's 
the incoming electron is no longer able to 'ionize' the bound electron (see e. g. T(4, .1)). At 
this point, both the product of the transmitted and reflected wave functions and the product 
of the 'bound' and 'free' electron wave functions on the 'dot' are maximal, giving a peak 
in the interaction between the free and the bound electron and therefore in the screening. 
This peak disappears as p increases. The peaked structure of the transmission versus k, 
which represents the peaks of T as function of the gate voltage on the 'dot', is very different 
from that found in the naive Coulomb blockade picture. ]7[ In the latter case, the distances 
between peaks would be equal to U\ We conclude that this naive picture fails for our exactly 
solved example. 

We now give more technical details. We consider only the singlet scattering wave func- 



tions, which are spatially symmetric, ty(xi,x 2 ) = ^(x 2 ,xi). At total energy E, we split 
\1/ into ^ = + ^s? where \l/o is the 'incoming' solution of TIq, with the same energy 
E + = (E + irj) (with rj -> + ), (H - E + )^ = 0. It then follows that 

*s(x 1} x 2 ) = UG E (x 1 ,x 2 ; 0, 0)tf (0, 0), (4) 

where G E is the two-particle Green's function of the interacting Hamiltonian, obeying 

(H — E + )G E (xi,x 2 ] x[, x 2 ) = —S(xi — x' 1 )5(x 2 — x' 2 ). (5) 

For the model Hamiltonian given by (p]) one can express Ge in terms of the two-particle 
Green's function of the non-interacting system, G° E : 

(Ho — E + )G° E (xi, x 2 ; x[,x 2 ) = —5(x\ — x' l )S(x 2 — x' 2 ). (6) 

Combining Eqs. (]5|) and @ yields 

G E (x 1 ,x 2 ;x' 1 ,x' 2 ) = G E (x 1 ,x 2 ;x' 1 ,x 2 ) + UG° E (0,0;x' 1 ,x 2 )G E (x 1 ,x 2 ;0,0), (7) 

and hence 

G E {x 1 ,x 2 ; 0, 0) = G%( Xl , x 2 - 0, 0)/[l - UG° E {0, 0; 0, 0)] = F E G%{x u x 2] 0, 0)/U. (8) 

This immediately yields G E {x\, x 2 ; x[, x' 2 ) and 

*(ar x , x 2 ) = %( Xl , x 2 ) + F E G%{x 1} x 2] 0, 0)* (0, 0). (9) 

The right hand side of (|9|) is determined solely by the eigenstates of the non-interacting 
Hamiltonian H,q. U appears only in the form Fe, which turns out to be quite important in 
determining the transmission characteristics. 

Instead of solving directly for the two-electron non-interacting Green's function G° E , we 
present this function in terms of the single-particle Green's function, g e (x,x'), of H. sp , 

(H sp (x) - e + )g t (x, x') = -6{x - x'), (10) 
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where e + = e + it], and rj — > + . g e (x,x') can also be written in terms of its spectral 
decomposition, 

9((xy) = E «£)«M 

where TC sp (pn = £n4>n- Writing the non-interacting singlet spatial wave functions as 

% m { Xl ,X 2 ) = (Mxi)</>m(X2) +<f>n (x 2 )(f) m (x 1 )) /2^ 1+S "^ 2 , (12) 

with energy E(n, m) = e n + e m , one can show that the spectral decomposition of the singlet 
non-interacting two-particle Green's function is 

Using the two spectral representations, we obtain 

1 r°° 

GP E {xi, x 2 ; x[, x' 2 ) = — / de g E - e (x\, ^i)%e(^2, 4) = 

7T J-oo 

i f°° 

— de g E _ e {x 1 ,x' 1 )g t {x 2 ,x' 2 ), (14) 

Z7T J-oo 

where the last equality follows from the Kramers-Kronig relations. 

We next calculate the quantum average of the current density operator at x = Xq, in the 
exact singlet state Noting that ^>{x\,x 2 ) is symmetric in x\, x 2 , this average is 

2eh r d 

j(x ) = S / dxidx 2 5(xi - xq)^*(xi,x 2 )- — ^/(xi,x 2 ), (15) 

m J dxi 

where \I/ is given by Eq. @. The explicit calculation of j now requires only integrals 
involving the non-interacting functions ^o(xi,x 2 ) and g € (x,0). In what follows, we shall 
assume that \l/ = ^o 9 > as given by Eq. (|12|), and that the total energy is given by E = e p +e q . 
The calculation is then facilitated using identities such as 

dx4," r (x)g,_(x,0) = ^(0)/(e+ - £„); 

e 2 — ei + ir) 

J deg E _ € (x, 0)/ (e + - e p ) = -2mg E _ ep (x, 0), (16) 



where the last equation represents the Kramers-Kronig relations. 

To proceed, we need to specify H sp . As the simplest possible example, we choose a simple 
(^-function attractive potential, 

which has one bound state, 4>b — \f^e~ K ^ with the inverse localization length k = mV/h 2 , 
and with eigenenergy — e& = —h 2 K 2 / (2m), and "band" scattering wave functions <p p = (e vpx + 
r p e w ^)/y/L, with the reflection and transmission amplitudes 

r p = in/(p-iK), t p = p/(p-iK), (18) 



and with eigenenergy e p = h 2 p 2 /(2m). For this simple Hamiltonian one has [14 



e (x,x') =g e (x',x) = -im(e lk ^- x '\ + r kc e ik ^ + ^) /(fi%), (19) 



where the wave vector k e is defined by k e = y/2me + /h, with Qk t > 0. 

There are two physical situations which are of interest for the non-interacting wave 
function ^q 9 . The first corresponds to two propagating electrons, impinging from the left, 
when both p = p\ and q = pi represent "band" states. This yields the simple Eq. (g). We 
devote the rest of this paper to the second, more interesting, choice for ^q 9 , when one electron 
is propagating, with p representing its wave vector, and the other is captured in the bound 
state "6", i. e. q = in. Now the total energy is E(p, b) = e p — e b = Ti 2 (p 2 — K 2 )/(2m) > —e b . In 
this case, the terms coming from the interaction are of the same order as the non-interacting 



ones, as now the amplitude for the two electrons to be together at x = is of order tJk/L. 
A long calculation, using Eqs. (pT3[), (|T4]), ([16]) and (|19D, now yields Eq. (|3|), and the effective 
transmission is given by 

T(U, p) = \t p \ 2 - 2mR(F E r p t p )/Ti 2 . (20) 

The second term here will have a "resonance" (yielding the large-K peak in Fig. 1) when 
the real part of the denominator of Fe (cf. Eq. @) will vanish, i. e. when 



U^E)- 1 = ^(Fe 1 ) = U- 1 - (0, 0; 0, 0) (21) 

vanishes. In some sense, U e s{E) represents the renormalized interaction. The details of the 
transmission thus require explicit expresions for G E (0, 0; 0, 0). Introducing an upper cutoff 
W = h 2 u 2 /2m on the "band" states, and using Eq. (0), we find (for real p > 0) 

1 K 2 , 1 IWs 

■— \ny c + -5- — -(1 arctanW — 

47T K z + p z K 7T V e & 



KG° s (0,0;0,0) = -2 
ft 



+ -s In — f= ; - In 



/t 2 +p 2 27rV _ + e b \l-y c yi\ 

2 



(22) 



with y c = (2yJW(W-E) + 2W - E) /E 2 , yi = (p + k) 2 /{ P - k) 2 , and 



9fG| 0, 0; 0, = --j. + Q(E) (£— ) - T — -J , 23 

where Q(E) is the Heavyside function. 

When p = 0, then G^(0, 0; 0, 0) is real and the equation dtG E (0, 0; 0, 0) = 1/U (equivalent 
to the resonance C/ e g {E)^ 1 = 0) is identical to the equation for the transition from a doubly 



bound ground state to a ground state with one 'free' electron. |jT2| In our special case, 
$IG E (0, 0; 0, 0) starts at a negative value at k — 0, changes to positive values around k ~ 
0.3, and approaches the asymptotic value 1/2 at large k. Thus, the above equation has 
a solution only for U > 2, when such a transition occurs. For small k, 3tG E (0, 0; 0, 0), 
remains negative, U e g(E) remains repulsive and there is no "resonance"; the transmission 
increases monotonically with U, as shown in Fig. @(a). For large k (and p > 0), U e g(E) 
becomes negative, reflecting an attractive effective interaction! An explicit calculation now 
yields a peak in T(U,p) just before U e g(E) changes sign, see Fig. @(b). These results are 
understandable qualitatively: at small U, the presence of the "bound" electron on the "dot" 
weakens the attractive potential V, and thus causes an increase in the transmission. Indeed, 
a Hartree-Fock-like approximation, in which we calculate the average of the interaction term 
with the symmetrized wave function \&Q b , yields V — ► V—2Uk, and thus k — > K(l — 2mU/h 2 ). 
Using this renormalized value of k in the bare transmission \t p \ 2 = p 2 /(p 2 + k 2 ), reproduces 
the result ( ffi^ ) in the small U limit. Thus, the interaction U compensates the (single- 
particle) attraction, and causes an increase in T(U,p). |L5| When C/ e g(£ l ) changes sign and 
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becomes negative, the transmission has a pronounced peak. A similar dramatic behavior of 
T(U,p) is manifested by its dependence on the depth of the single-particle bound state, V, 
or - equivalently - on k, as shown in Fig. |I[ 

The above results indicate that the more interesting behavior of the transmission occurs 
when the single-particle attraction is high enough, k > u (that is, when the bound state 
energy |e&| is larger than the band width, W). In this situation, the total energy of the two 
electrons is negative [p 2 < u 2 < k 2 }, and QG° E is given only by the first term in Eq. (]23|). 
Then, the transmission takes an especially simple form 

T(U,p) = \t p \ 2 , t p = t p (l - F E r p ) = + M 2 )/(^ + i-jr^). (24) 

Hence, when 1/C4ff approaches 0, t p — > —r p . At large values of k the reflection is close to 
unity. This explains the heights of the peaks in Figs. 0(b) and [I]. 

Returning to Eq. (|22|), is is interesting to note that the sign change in 31G E originates 
from the last two terms, which vanish when k = 0. In the absence of the bound state, U e s 
always remains negative. In the usual theory of superconductivity, the remaining first term 
in Eq. (^2p renormalizes U into a weaker repulsion, crucial for the superconducting state. 
All the interesting phenomena found here result from the additional effects of the bound 
state. 

Finally, a comment about spins: starting with two general spins o~\ and o"2, one can 
always split the wave function into a combination of a singlet and a triplet. Without the 
interaction, the scattered electrons will have the same spins. Given the above results, the 
interaction will change the relative weights of the singlet and the triplet, and thus may cause 
a spin flip. 
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FIG. 1. The transmission T(U,p) as function of k (in units of the cutoff w), for U = 1, 2, 4 
(in units of fi 2 /{2m)), solid line, and T(0,p) = |t p | 2 , dashed line, p is in units of u>. 
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FIG. 2. The transmission T(U,p) as function of the 'bare' interaction U (in units of ft 2 /{2m)). 
The solid line is T(U,p), the dashed line shows \t p \ 2 . (a) $tG° E < (k = 0.2a;, p = 0.04a;). (b) 
$IG E > (k = uj, p = 0.1a;). The thick point signifies the change of sign of U e g. 
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